celar

*Figure 1

use "data_for_figures.dta", clear

twoway (kdensity dds_mom if treat_all_round==0 & round==2, lwidth(thick) lcolor(black) lpattern(solid)) (kdensity dds_mom if treat_all_round==1 & round==2, lwidth(thick) lcolor(gs7) lpattern(dash)), title("{bf: A}", size(large) ring(1) pos(11) span) xtitle("Mothers' dietary diversity", size(medlarge)) ytitle("Fraction", size(medlarge)) ylabel(0(0.1)0.5, labsize(medium) nogrid) xlabel(0(1)8, labsize(medium) nogrid) legend(off) graphregion(margin(zero)) saving(Baseline, replace) 
twoway (kdensity dds_mom if treat_all_round==0 & round==3, lwidth(thick) lcolor(black) lpattern(solid)) (kdensity dds_mom if treat_all_round==1 & round==3, lwidth(thick) lcolor(gs7) lpattern(dash)), title("{bf: B}", size(large) ring(1) pos(11) span) xtitle("Mothers' dietary diversity", size(medlarge)) ytitle("Fraction", size(medlarge)) ylabel(0(0.1)0.5, labsize(medium) nogrid) xlabel(0(1)8, labsize(medium) nogrid) legend(label(1 "Control") label(2 "Treatment") pos(3) ring(0) cols(1) size(medsmall) ) graphregion(margin(zero)) saving(Post-treatment, replace) 

graph combine Baseline.gph Post-treatment.gph, cols(2) rows(1)
graph display, ysize(2) xsize(4) scale(1.3)
graph export "dds_mom_fatigue.eps", replace